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Dynamic and thermodynamic magnetization experiments on heavy-ion irradiated single crystalline 
Bi2Sr2CaCu208+5 are correlated in order to clarify the nature of the mixed state phase diagram. 
It is shown that whereas the entropy contribution to the free energy in the London regime plays a 
minor role in unirradiated crystals and irradiated crystals at fields close to or above the matching 
field B^, it becomes very important at low fields in irradiated crystals with high B^. The direct 
determination of the entropy contribution to the free energy from the reversible magnetization 
allows one to determine not only the correct values of the pinning energy, but also to extract quite 
detailed information on pancake vortex alignment. The characteristic field Hint ~ gB^ at which 
intervortex repulsion begins to determine the vortex arrangement and the reversible magnetization 
is shown to coincide with a sharp increase in the irreversibility field Hirr{T) and with the recoupling 
transition found in Josephson Plasma Resonance. Above Hi„t, the repulsive interaction between 
vortices cause both the vortex mobility to decrease and pancake alignement to increase. At higher 
fields > \B,f, 2> -Bel, free vortices outnumber those that are trapped on a columnar defect. This 
causes the decrease of c-axis resistivity and a second crossover of the irreversibility field, to a regime 
where it is determined by plastic creep. 
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I. INTRODUCTION 

The remarkable influence of heavy ion-irradiation in- 
duced amorphous columnar defects on the mixed state 
phase-diagram and transport properties of layered su- 
perconductors, and the interplay between the pinning 
of vortices and couphng betweea. the layers, has pro- 
voked a great deal of interest .Ert3 The effect of the ion 
tracks is twofold. First, the destruction of supercon- 
ductivity within the columnar defect core, the radius 
of which is comparable to the superconducting coher- 
ence length ^, leadsj-^o a large reduction of the mixed 
state free energyfij t3 arising from the pinning of the 
two-dijnensional pancake vortices in the superconducting 
layergl3E3 by the columns. The reduction in free energy 
due to pinning is readily measurable because it entails 
a significant reduction of the reversible magnetization 
Mrev at fields well below the matching field i?^ = '^arid, 
and a nonmonotonicity of dMrev/dlnB as function of B 
(<I>o = h/2e is the flux quantum and Jirf, is, the density 
of column intersections with the layers). L3t3 The coinci- 
dence of the radius of the columnar defect cor|e-|With that 
of the vortex core implies a large pinning forccEil and coru 
committantly large values of the critical current density,El 
as well as a large increase of the field Hirr [T) below which 



irreversible magnetic response is measuredJa Second, the 
linear geometry of the ion tracks results in the effective 
alignment of the pancake vortices along the tracks, and 
the re-establishment of interlayer superconducting phase 
coherence at temperatures and fields where the unirradi- 
ated material behavca.a.a_a .stack of nearly decoupled su- 
perconducting layers .qQ'Ej Ej The pancake vortex align- 
ment is manifest in the angular dependence of the tfAns- 
port properties of the irradiated superconductor,Q'El as 
well as in measurements of the c-axis condufjtivitjy (per- 
pendicular to the superconducting layers)BEjc3'LZl and of 
the field at which the Josephsoti. xJasma resonance oc- 
curs in the range 24 - 45 GHz.EJ'E3l30 Measurements 
of the JoscnhsQii plasma resonance and of the c-axis 
conductivity,Ejc3 which have the advantage of directly 
probing the Josephson critical current and the cosine of 
the average phase difference between layers, show the ex- 
istence, in the heavy-ion irradiated layered superconduc- 
tor BiaSraCaCuzOg+s, of a regime Bi{T) < B < Bh{T) 
in the vortex liquid state in which interlayer supercon- 
ducting phase coherence increases as the magnetic field 
is increased {i.e. a "recoupling"). This behavior is quali- 
tatively opposite to that found in the unirradiated mate- 
rial. At fields greater than Bh, interlayer phase coherence 
again decreases as field is increased. 



evcrsiblc 
It is the 



The peculiar field dependence of the c-axis electromag- 
netic response, as well as that of the reversible - and 
even of tJae irreversible part of the magnetization (below 
Hirr{T))l3 have been interpreted in terms of the com- 
petition between the Josephson coupling energy, which 
tends to align pancake vortices, and the entropy con- 
tribution to the free energy arising from a random dis- 
tribution of pancakes over columnar defects. lj The im- 
portance of intralayer pancake repulsion is often down- 
played, even though it has been shown to Ipe, important 
in determinina— |i|iip pancake arrangement 
magnetizationJU't^l and transport properties 
purpose of this paper to investigate the importance of 
the contributions of intra- and interlayer pancake vortex 
interaction energy, entropy, and pinning energy to the 
Gibbs free energy G, and the impact of each on the phase 
diagram of heavy-ion irradiated Bi2Sr2CaCu208+5- To 
this effect, we have performed reversible magnetization 
as well as AC screening measurements on heavy-ion irra- 
diated Bi2Sr2CaCu208+5 with widely varying matching 
fields (section ||). From the analysis presented below 
(section [II A| ) , it turns out that in the London regime, 
the entropy contribution to the reversible magnetization 
is rather insignificant with respect to the total magneti- 
zation, or its modification due to vortex pinning on the 
columnar defects. Only in irradiated crystals with suf- 
ficiently high i?0 is it essential to take a configurational 
entropy contribution into account when describing the 
magneti zation, and then, only at low fields B <C B^ 



(section |III A2| ). Asj-ajjaDnsequence, previous estimates 
of pinning energiesOlljtj determined from the reversible 
magnetization should be revised towards (sometimes sig- 
nificantly) lower values. Next, it is shown how informa- 
tion on pancake alignment can be obtained from the re- 
versible magnetization (section III A 3). The data allow 
for the extraction of the field Hint at which vortex repul- 
sive interactions begin to determine the pancake vortex 
arrangement over the columnar defects. This field corre- 
lates very well with sharp features in the irreversibility 
line, and, in some cases, with the "recoupling" observed 
in the JPR measurements and c-axis conductivity ( IIIBJ ). 



II. EXPERIMENTAL DETAILS 

A. Sample preparation 

For this study, we have used Bi2Sr2CaCu208+5 sin- 
gle crystals grown at the University of Amsterdam by 
the travelling solvent fioating zone techniqueS A num- 
ber of crystals were postannealed at 800°C in air and 
had Tc = 90.0 K. Other crystals were retained as-grown 
— this yields samples that are lightly overdoped in oxy- 
gen, and that have Tc « 83 K. Another batch of crys- 
tals, groW|tL.at the University of Tokyo using the same 
technique J2j had Tc = 88.8 K and was used for the mea- 
surements of the irreversibility line for different match- 



ing fields. After growth and annealing, large crystals 
were selected for uniformity and absence of macrodefects, 
and subsequently cut using a wire saw, so as to produce 
small squares of typical size 800 x 800 x 20 /im"^, suit- 
able for magnetometry experiments. A number of larger 
pieces (1x2 mm^ x 20 /im) of the Amsterdam crystals 
were kept for reversible magnetization measurements us- 
ing a commercial superconducting quantum interference 
device (SQUID) magnetometer, as well as ac transmit- 
tivity measurements. In this manner, a sizeable number 
of crystals with similar characteristics was obtained. 

The samples were irradiated with 5.8 GeV Pb ions at 
the Grand Accelerateur National d'lons Lourds (GANIL) 
at Caen, France. Such an irradiation produces continu- 
ous amorphous tracks of radius cq ~ 3.5 nm which tra- 
verse the samples along their entire thickness. During 
the irradiation, the ion beam traverses a 1 /im thick Ti 
film placed in front of the sample; the secondary electron 
emission from this film is used to continuously monitor 
the ion flux during irradiation. At the beginning of each 
run, the ion flux is calibrated using a Faraday cup placed 
between the Ti-film and the sample. This procedure al- 
lows one to accurately control the total ion fiuence for 
each sample. Moreover, the ion beam is continuously 
swept across the target using asynchronous vertical (3 
Hz) and horizontal (1 kHz) ac drive fields so as to ex- 
pose the entire target area (2x3 cm^) homogeneously. 
Different crystals from each source were irradiated with 
widely different ion fluences, in order to produce samples 
with dose equivalent matching fields 0.02 T< B^ < 4 T. 
There was a slight reduction of the samples' Tc after the 
irradiation. Notably, overdoped single crystals had a Tc 
of 81 K and 78.7 K after irradiation to doses correspond- 
ing to B^ = 2 and 4 T respectively. The crystals from 
the University of Tokyo underwent a reduction of Tc of 
up to 4 K (for S^ = 4 T). 



B. ac shielding 

The ac transmittivitytJ was measured using the Lo- 
cal Hall Probe Magnetometer, in the same way as de- 
scribed in Ref. ^. The sample was placed on top of a 
miniature InSb Hall sensor; both were surrounded by an 
AC drive coil that could produce fields hac of up to 30 
G (at frequencies 0.5 Hz < / < 2 kHz). The ac field 
was always directed parallel to the sample c-axis. This 
corresponds to the application to the sample of an az- 
imuthal electric field E^ with gradient r^'^d{rE^)/dr — 
2'Khacf {r is the radial coordinate). The ac Hall volt- 
age Vac was measured as function of temperature and 
applied field Ha- From Vac, the first and third har- 
monic components of the transmittivity are determined 

as TH{f,T) = {Vacif,T) - Vacif,T<^Tc))/iVac{f,T-> 

Tc) - Vac{f,T^Tcj) and Ths = Vac{3f,T)/{Vac{f,T:^ 
Tc) — Vacif, T<^Tcj) respectively. In order to determine 
the irreversibility field Hirr{T, /) (or irreversibility tem- 



perature Tirr{Ha, f)), the transmittivity was measured 
using an ac field of amplitude hac = 1 G and frequency 
/ = 7.75 Hz. Tirr{Ha, f) was taken as the temperature at 
which the Tf/3 signal first becomes distinguishable from 
the background noise when cooling. Given the size of our 
samples, this corresponds to probing whether the E{j)- 
characteristic deviates from Ohm's law at an electric field 
value of 0.6 nVcm~^. In practice, the measurement of 
Tirr is often limited by instrumental resolution, which 
can be estimated by comparing the signal at full screen- 
ing to the noise level; it corresponds, for a sample of area 
800 X 800 jiTD? , to a circulating shielding current density 
of 4x 10^ Am-2. 



C. Reversible magnetization 

Measurements of the reversible magnetization were 
carried out using a commercial SQUID magnetometer 
at the University of Leiden. Care was taken to avoid 
the numerous artefacts previously observed in such mea- 
surements. Noteably,|-w'e have found that the reported 
"clockwise" hysteresisE3 is an artefact related to the use 
of a too restricted length over which the sample is drawn 
through the SQUID gradiometer coils when the magnetic 
moment is small. A scanlength of 6 cm produced satis- 
factory results, with a reproducible reversible magnetic 
moment at fields greater than Hirr{T). However, it is 
well-known that, as a consequence of the inhomogeneity 
of the field produced by the superconducting solenoid, 
measuring the magnetization with such a scanlength at 
fields close too but below Hirr corresponds to cycling the 
sample through minor magnetic hysteresis loops. The 
magnetic moment then results from the superposition of 
several ill-controlled shielding current loops and can take 
on a value that is very different from that expected from 
the ordinary full critical state. It thus proved very dif- 
ficult to accurately measure the magnetic moment both 
above and below Hirr{T); for the present paper, we chose 
conditions such as to obtain the most reliable results 
above Hirr (T) , and to rely on the ac transmittivity mea- 
surements to obtain the screening current below Hirr(T). 



III. RESULTS AND DATA ANALYSIS 

A. Reversible magnetization 

Figure P shows the reversible magnetization for an op- 
timally [doped Bi2Sr2CaCu2 08+5 single crystal {Tc — 
88.8 K),E3 irradiated with 10^^ Pb ions cm~^, which cor- 
responds to a matching field B^ = 2 T. This valuo-af i?^ 
is comparable to that used in previous studiesJ13il3 and 
the reversible magnetization shows all the features dis- 
cussed there. Recapitulating, at low fields Hirr < Ha -^ 
B^, the absolute value of the magnetization |Mret;| de- 
creases more or less proportionally to \n{l/ Ha), in 




^H(T) 

FIG. 1. (a) The reversible magnetization of a heavy-ion 
irradiated optimally doped Bi2Sr2CaCu2 08+a single crystal 
{B^ = 2T;Tc = 88.8 K), for temperatures T = 70, 74, and 
78 K. The drawn lines represent fits to Eq. (H) with parameter 
values Al(0) = 240 nm (A(0) = 180 nm), dBc2/dT = -1.05 
T/K, and Uqs ~ l.bSos r; 1300(1 - T/Tc) K. Arrows indicate 
the field below which the ac transmittivity has a non-zero 
third harmonic (Hirr) and the interaction field Hi„t above 
which vortex interactions infiuence the vortex arrangement; 
(b) (c&os/fcs = 0.23 Km/A) times the difference between 
the measured magnetization and the theoretical magnetiza- 
tion in the absence of columnar defect pinning, M^eu (cor- 
rected for the reduction in superffuid density expressed by 
A"^(T) = (1 - 27rcond)A"^). In the limit B < B4,, the plot- 
ted data correspond exactly to the sum of pinning energy 
and entropy per pancake vortex. Horizontal dotted lines cor- 
respond to the bare pinning energy L'^os ~ t^o(0)s(l — T/Tc)^ , 
with C/o(0)s — 4000 K. This corresponds to the parameter 
value obtained in Ref. h_q, corrected for the change A ^ A. 
The drawn lines show the logarithmic extrapolation of the 
low field data which permit the determination of the length 
La over which pancakes belonging the same vortex line are 
aligned on the same column, and the fac tor d escribing the 
fraction of available columns (see section p!II A3 ) . For T = 70 
K, the low field data obey Uqs + T\n{f3B^/B), i.e. the typi- 
cal stacklength La is close to the interlayer spacing; for larger 
temperatures the prefactor to the logarithm is, unexpectedly, 
somewhat less. The interaction field Hi„t of Fig. hl(a) is deter- 
mined as that where the data first deviate from a logarithmic 
Ha-dependence. 

agreement with the London model which has Mrev ~ 
-{eo/2<^o)HvBc2/eB) = M°^,. Here £0 = 

<i>o/47r^oA^(r), A(T) is the penetration depth, Bc2 is the 
upper critical field, ^q = Air x lO^^Hm^^, i] ^ 1 and 
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FIG. 2. Logarithmic slopes dMrev/d In B for the same crys- 
tal as in Fig. nl, determined in the limits B > B^ = 2 T 
and B <^ B^ respectively. The drawn lines indicate the 
magnitude of the theoretical entropy contribution in the 2D 
hmit, fcsT/$os, and the prediction of the London model, 
M°e„ =£o/2$o. 

e = 2.718 . . .. Throughout, the small magnitude of Mrev 
allows us to take B « HoHa- The logarithmic de- 
crease of lAfreul terminates at the field labelled Hint, 
whence |Mret)| increases. When B ^ i?^, |Afret)| de- 
creases once again proportionally to ln(l/i3), although 
the slope dMrev/dlnB is markedly higher than the slope 
at low fields. The logarithmic slopes in either limit 
are illustrated in Fig. g, which shows dM^ev/dhiB as 
function of temperature for the same crystal: for fields 
B ^ Bfp the logarithmic slope lies markedly below that 
measured for B > B^; also, it goes to zero at a tem- 
perature Ti that is lower than the temperature T2 at 
which dMrev{B > B^)/dB goes to zero. Both T2 and 
dMrev{B > B^)/dB are very similar to the correspond- 
ing values measured in the unirradiated crystaLlij In 
both cases, unirradiated crystals and irradiated crys- 
tals in_tlip high field limit, one can use the London 
modeOcHl to extract the value of the penetration depth 
from d^Mrev{B > B^)/dTd\nB = -eo(0)/2$oTc, an 
exercise which yields A(0) = 180 ± 30 nm. The corre- 
sponding value of the London penetration depth Al (0) ~ 
1.35A(0) w 240 nm is in good agreement with|4iterature 
values for optimally doped Bi2Sr2CaCu208+5E^ The er- 
ror (of less than 20 %)cj arises because of two reasons. 
First, the presence of the normal cores of the colum- 
nar defects decreases the supercurrent density around 
each vortex and the strength of intervortex repulsion. 
This effect can be taken into account by a "renormal- 
ized" value A~^(T) = (1 — 2TTCQnd)X~^ and correspond- 
ing io = $g/47r/ioA2 and M,% = eo/2<^oHvBc2/eB)E 
More important, the presence of thermal fluctuations was 
predicted to lead to a substantial entropy caatribution to 
the free energy, reflected by an extra term:c3 



dTdlnB 
is half the 



(1) 



dHdre^ _ dHl^ _ kB_ 

c-axis parameter of 



dTdlnB 

Here, s ~ 1.5 nm 
the Bi2Sr2CaCu208-i-5 material {i.e. the distance be- 
tween Cu02-bilayers). From the value ks/^os = 4.4 



Am~^K~^, we-see that the (constant) contribution of the 
entropy termcj can amount to 20 % of the temperature 
derivative of dMrev/dlnB, at most. 

The nonmonotonic field dependence of the reversible 
magnetization has a straightforward explanation in terms 
of the lowering of the wjttex free energy by pinning onto 
the columnar defects. t30 The magnetization Mrev — 
-dG/dB = ~^Q^dG/dnv = -^q V: with Uy the vor- 
tex density, corresponds to the vortex chemical potential 
/i and therefore directly measures the energy needed to 
add a vortex to the system when the field is increased. At 
low fields many columnar defects are available, so that 
each and every new vortex can gain a maximum energy 
per unit length Uq by becoming trapped in the defect 
potential; correspondingly Mrev ~ Mrev + C^o/^o- The 
pinning of a vortex on an insulating columnar defect en- 
tails the redistribution of the supercurrent towards the 
vortex periphery, lower maximum current on the vor- 
tex core boundary, and a smaller intervortex repulsion. 
More vortices can therefore enter the sample at fixed 
magnetic field, hence the absolute value of the magne- 
tization is smaller. However, it . i s. a i|4U expected that 
Mrev oc ln(l/B) in this regimcNoH At high fields 
fioHa > B^ nearly all columns are occupied, new vor- 
tices cannot be trapped on a defect, hence Mrev ~ M^ev^ 
close to the magnetization of an unirradiated sample in 
the London regime.cJ In the intermediate field region, 
Mo-ffjnt < ^J■oHa < Bff, in Fig. |l|(a), new vortices can be- 
come either trapped or free. The proportion of trapped 
vortices decreases as field increases, leading to the in- 
crease of |Mret, |. 



1. Role of intervortex repulsion 

The magnetization in the intermediate regime has been 
estimated in Ref. |^ by considering that the gain in pin- 
ning energy for occupation of a particular column should 
be larger than the loss in energy due to intervortex re- 
pulsion when a vortex is moved from its lattice site to 
the columnar track; assuming a Poisson-distribution of 
distances between tracks, this gives 



Mrev ~ M\'.^y 



Uo 



Uo B^ 
1^~B 



-{Uo/io){B^/B) 



(2) 



In spite of its simplicity, it is possible to obtain accept- 
able fits to the experimental data at Ha ^ -ffmt us- 
ing this expression (see Fig. ma) and Ref. B2h. Jji the 
Figure we have used dBe2/dT = -1.05 T/K,Ea and 
A(0) = 180 nm as determined above. However, for a 
faithful fit to the peak structure of |Afret)| one should 
assume Uq/eo to be in excess of unity; for example, at 
70 K, Uq/sq = 1.5 in Fig. ||. Such large values of Uq 
were also found in Refs. O and p2; they are incompati- 
ble with the theoretical expectation for electromagnetic 



pinning of vortices,E3ea'c2l Uq w eq ln(co/\/20' since they 
would imply that the column radius cq exceeds the co- 
herence length ^ = 2.2(1 - T/Tc)-^/'^ nm by a factor 5 - 
6. In Ref. |l8|, it was found that pinning of the vortex core 
dominates over electromagnetic pinning in heavy-ion ir- 
radiated Bi2Sr2CaCu208+5 at all temperatures. Then, 
Uo « £o(co/2^)2, yielding cq « 2.3^(0). This value is 
still somewhat high, which shows that the actual pin- 
ning mechanism is more complicated than simple core 
pinning. However, note that even so the magnetization 
data at fields Ha <C Hint lie much above the theoretical 
curve. 

We thus find that a description that only takes the 
decrease of the internal energy due to the localization 
of the vortices on the columnar defects into account is 
insufficient to describe the large decrease of the low-field 
magnetization at large S^. The fact that the low-field 
logarithmic slope dM/dlnB is much smaller than that 
predicted by Eq. (Q) indicates the presence of an extra 
contribution to the free energy that has a logarithmic 
dependence ^ \n{l/B). 



2. Role of entropy 

In this section we examine the importance of a possi- 
ble entropy contribution —TS to the free ene rg j\ r , i Su ch a 
contribution may arise from vortex positpnalr'y4c3 and 
order parameter amplitude fluctuations O'tJO The en- 
tropy contribution to the free energy is believed to be im- 
portant in layered superconductors because of the "cross- 
ing point" behavicp-o b i- icij | i . rq d in the reversible magnetiza- 
tion at T* and T^ MtjEaH Supposing that a description 
in terms of the London model is appropriate, and hence, 
that vortex positional fluctuations give the dominant con- 
tribution to S, it was derived in Ref. BS that the crossing 
point arises in the unirradiated material when the field 
dependence of the two terms contributing to Mrev , 
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cancel at T* = Ta/{1 + 2kBTc/eo{0)s). ^^Fi^o is a field 
determined by the elemental phase area;L2lc3 the fraction 
Bq/B corresponds to the total number of possible states 
( "configurEitipns" ) per vortex. In two dimensions (2D), 
Bq ^ Bc2-'^ In heavy- ion irradiated layered supercon- 
ductors, Eq. (0) should be modified: the full result can 
be found in Ref. Bil but the limiting forms for Mrev are: 



Mrev « Ml!^^ 



Mrev « Ml!^^ 
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{B<^B^) (5) 
(B»S^;[/o^O). (6) 



These are again valid in the 2D limit where the posi- 
tions of pancake vortices in neighboring layers are uncor- 
related. The logarithmic field dependence of the second 
term of Eq. (|^) now depends on the number (~ B^/B) of 
columnar defect sites available to each pancake vortex. 

In order to estimate the importance of the entropy con- 
tribution to the experimentally measured magnetization, 
we return to the temperature dependence of the logarith- 
mic slopes dMrev/dlnB in Fig. 0. Note that, according 
to Eqs. (|^-|^), dMrev/d\nB is equal to the difference be- 
tween two large terms: 



dMr 



dAlL 



kuT 



dlnB dlnB $os 



(7) 



However, these terms have very different temperature 
derivatives, equal to — eo(0)/2$oT'c and —ks/^os respec- 
tively. The ratio between the temperature derivatives 
is equal to 2Gi^'^, where Gi is the 2D Ginzburg num- 
ber that determines the regime of reduced temperature 
t{T) = {Tc — T)/Tc over which the entropy contribu- 
tion to the free energy is dominant. For optimally doped 
BiaSraCaCuaOg+A-, 2kBTe/eo{0)s « 0.19, i.e. Gi « 0.01; 
for our overdoped material, Gi « 0.003. The value 0.19 
means that the error in the A(0)-value deduced from the 
reversible magnetization, arising from an entropy contri- 
bution to the free energy, cannot exceed 19% (11% for 
the overdoped crystals). Inserting the experimental A- 
value back into Af"g^, and comparing the result to the 
full measured magnetization, we|fipd that for B ^ B^ as 
well as for unirradiated crystals,CJ the entropy contribu- 
tion to dMrev/dlnB does not exceed 0.1 ksT/^os (see 
Fig. 0). This excludes an interpretation of the crossing 
point in terms of quasi~2D vortex positional fluctuations 
only: for example, Eq. (||) would lead to T* — 75 K. 
More generally, vortex translational degrees of freedom 
do not lead to a major modiflcation of the magnetization 
of unirradiated crystals (or irradiated crystals at fields 
above i?^) in the London regime. The entropy contribu- 
tion to the free energy only starts to be important in the 
critical fluctuation regime. 

With respect to this, we note that the analysis of 
Ref. ^ has demonstrated that the crossing point in 
unirradiated BiaSraCaCuaOs+a crystals indeed lies in 
the regime of (2D) critical fluctuations. An adequate 
description, vahd for B > \t{dBc2l dt)T=Ta^ was ob- 
tained using in Ref. Eq by expanding the order parame- 
ter in terms of lowest Landau level (LLL)-eigenfunctions. 
The expression for Mrev derived in Ref. |4^ (which also 
predicts dMrev/dlnB -^ £o/2*o for t > t{T*)) well 
describes the experimental data for both unirradiated 
BiaSraCaCuaOs+^a ancL—irradiated BiaSraCaCuaOs-i-ij 
above the matching fieldJlj In Fig. 0, the B 3> -B^-data 
lie outside the London limit for T > SOisT, and within the 
LLL-regime for T > 84 K. An interpretation of the field- 
independent magnetization at Tj* in terms of vortex posi- 
tional fluctuations (order parameter phase fluctuations) 
only is thus inappropriate. 



The situation is totally different for B <C B^f,, at which 
one has a significantly larger value TS > 0.3kBT/^oS- 
This means that the introduction of many extra possible 
(columnar defect) sites within a vortex lattice unit cell 
leads to a large (configurational) entropy contribution. 
Moreover, a description in terms of the London model 
will turn out to be adequate in this field regime. 



3. Interlayer coupling: dependence on matching field 

The small entropy contribution to the reversible mag- 
netization for B > B^ suggests that even in the vortex 
liquid, significant positional correlations remain between 
vortex pancakes in adjacent layers. In that case, S is 
reduced with respect to its value in the 2D limit: the 
alignment of pancakes in different layers means that the 
relevant entities are stacks (of length La) comprising sev- 
eral pancakes. The entropy contribution in the regime 
B <^ B^ where all pancakes are trapped by a column 

can be estimated as TS = ksTlnW = fcsTlniVi'^/^", 
with Nc = PB^/B the number of columns available to 
each stack of pancakes, Lg the length of the sample in 
the field direction, and La the length \\ c over which the 
pancakes are aligned on the same columnar track. The 
number 1 — /3 is the fraction of columns that is inac- 
cessible due to intervortex repulsion. In the quasi-2D 
limit. La = s, but in general La can be estimated from 
the balance between elastic energy and thermal energy: 
ii^o/B^La ^ ksT (ei « £0/7^ is the line tension and 7 
the anisotropy parameter). Hence the entropy contribu- 
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{B^ < B„eos/kBT). 



The crossover field Bcr = ^i^/ilsY delimits the (low 
field) regime where tilt deformations of the vortex lat- 
tice with wavevector ~ 47r(i?/$o)^^^ are more difficult 
to excite than shear deformations with similar wavevec- 
tor, from the high field regime where the inverse is true. 
The contribution TS adds a supplementary logarithmic 
dependence on B at low fields, in agreement with the dif- 
ference between the experimental dMrev/d\a.B observed 
at i? <C B^ and B > B^ respectively (see fig. 0(a)). 

In Ref . 03, it was claimed that the maximum in | Mrev \ 
can be nearly entirely explained by the reduction in TS, 
caused by pancake alignment ( "recoupling" ) and the con- 
commitant decrease in s/La- This claim can be tested 
by using Eq. (Q), which permits a direct estimate of 
the correlation length La- Starting from the values of 
dMrev/d\nB determined above, one would have 




FIG. 3. (a) Reversible magnetization of lightly overdoped 
Bi2Sr2CaCu208+5 at T/Tc = 0.89 , for matching fields 0.2 T 
< B,jy < 4 T . The critical temperature Tc ~ 83.1 K, except 
for 50 = 2 T (Tc = 81 K ) and 50 = 4 T {Tc = 78.7 K). 
The drawn lines through the data points for B^ = 0.2 and 0.5 
T represent fits to Eq. (0) with parameter values A(0) = 140 
nm and Uqs = 120 K; the dotted lines are guides to the eye, 
representing the magnetization in the irreversible regime, (b) 
^os/kB times the difference between the measured magne- 
tization and the theoretical magnetization in the absence of 
columnar defect pinning, M^ev , (corrected for the reduction in 
superfiuid density expressed by A~^(T) = (1 — 27vCQnd)\~^). 
In the limit B ^ B^, the plotted data correspond exactly 
to the sum of pinning energy and entropy per pancake vor- 
tex. The drawn lines show the bare pinning energy per pan- 
cake vortex Uqs (corresponding to the low field limit of the 
-B0 = 0.5 T -data), and, for B^ > 1.0 T, C/qs plus the config- 
urational entropy contribution TS ~ (fcsT/^os) ln{pB^/B) 
arising from the possibility of individual pancake vortices to 
occupy different column sites. The fraction /3 ~ 1/2 can be 
extracted from the intercept with the horizontal Uos-line, and 
is indicated by the bold arrows. The fields Hint where inter- 
vortex interactions begin to infiuence the magnetization can 
be determined as the field of first deviation from a logarith- 
mic field dependence and are indicated using the thin arrows. 
Dotted lines are guides to the eye, indicating data in the ir- 
reversible regime. 

La < s/0.3 « 3s for ^o-^a < 50 == 2 T, and La > 10s for 
fJ-oHa ^ B^, similar to the value La ~ 15s recently ob- 
tained by Morozov et al.JEH and in qualitative agreement 
with the assumption of Ref. O. 

A more stringent approach is to consider the de- 
pendence of Mrev on matching field, i.e. the depen- 
dence on the number of available states. Figure ^(a) 
shows the magnetization of a series of lightly overdoped 



BiaSrsCaCuaOg+s crystals (T^ w 83 K, A(0) « 140 
nni), irradiated to different ion fluences corresponding 
to matching fields ranging between 0.2 T and 4.0 T. 
The Figure shows that there is a large dependence of 
the low-field magnetization on the density of columns, 
which irrevocably indicates the importance of the con- 
figurational entropy contribution to Mrev in that regime. 
Namely, if S were zero, the low-field limits of M^ev would 
have been equal to M°g^ + [/q/^o, independent of B^ 
(see. Eqs. (@,||)). Furthermore, the data at B^ > 1 T 
confirm the extra logarithmic field dependence (|8|) for 
B <^ B^. For the crystal with B^ = 0.5 T, the low- and 
high field limits of dMrev/ din B are nearly equal, imply- 
ing that at this temperature, the entropy contribution 
to d Mrev /d In B lies wi thin the error bar of O.l/csT/^os 
found in Section III A 2 (see also Ref. n% Fig. 3). Hence, 
for this matching field and temperature, it follows from 
Eq. (^ that pancakes are aligned into stacks of length 
La ^ 10s over the entire considered field range. It means 
that for low matching fields it is allowed to ignore the en- 
tropy contribution. One can then confidently determine 
UqS from the difference between the rBG|asured magneti- 
zation and the low field limit of M^^^ .llaH Inserting the 
obtained value, Uqs !=^ 120 K (for ios « 160 K at 74 K) 
back into Eq. (||) yields good fits to the data for i?^ = 0.2 
and 0.5 T, as shown by the drawn lines in Fig. y. Such 
fits are not possible at higher B^ due to the increasing 
importance of the entropy contribution, which, even at 
B^ = 1 T already accounts for half the change of the 
low-field reversible magnetization with respect to that 
of an unirradiated crystal. The determination of a "pin- 
ning energy" from the difference between the low field 
data and the extrapolation to low fields of the data for 
B ^ i?0 can then introduce an error of more than a 
factor 2. 

The parameters La and /? can, in principle, be di- 
rectly evaluated by subtracting M^^^ from the measured 
magnetization. This procedure yields a quantity that at 
low fields is rigorously equal to the sum of the pinning 
energy per pancake vortex and temperature times en- 
tropy (Fig. 1(b)). The low-field fimit of the B^ = 0.5 
T data corresponds to the bare pinning energy Uqs — 
this value is indicated by the horizontal line in Fig. y(b). 
The results for the crystals with B^ > 1 T exceed Uqs 
at B <C B^ because of the entropy contribution to the 
free energy. From the logarithmic field derivative and 
the extrapolation to Uqs we find that, at T/Tc = 0.89, 
La{B <C Be/,) « s and /3 ~ i, independent of B^. Thus, 
the magnetization at low fields and B^ > 1 T is well 
described by 



Mrev = M°ev + ^ 
(f>0 






(10) 



as illustrated by the drawn lines through data in 
Fig. p|(b)). The first deviation of the experimental data 
from these lines as field is increased, at the "interaction 
field" Hint, indicates that both configurational entropy 



and the total pinning energy start to decrease due to the 
effect of vortex interactions. Entropy decreases because 
of the limitation of accessible column sites, which is the 
consequence of intraplane pancake repulsion; it can be ex- 
pressed by a diminishing /3; the translational invariance 
of the system along the column direction and the attrac- 
tive interaction between pancakes in adjacent layers will 
then lead to pancake vortex alignment and the supple- 
mentary decrease of S through the increase of La- The 
total pinning energy Up can only go down appreciably if 
free (or "interstitial") vortices, not trapped on a colum- 
nar track, start to appear, again as a result of intraplane 
pancake repulsion which prohibits certain vortices from 
finding a favorable column site. From Fig. Bl it is seen 
that near Hint, the pinning energy (obtained from the 
B^ = 0.5 T data) and the entropy contribution to the 
free energy are comparable in magnitude when B^ ^ 1 
T; hence both contributions should always be taken into 
account when describing the shape of the Mret)(B)-curve. 

We have attempted to carry out the same analysis for 
constant matching fields B^ = 0.5, 1, and 2 T, and dif- 
ferent T (see Fig. ||(b) for B^ = 2 T). In this case, the 
problem arises that the bare pinning energy cannot easily 
be established. It was therefore chosen to take the result 
of Ref. |l| (where B^ = 0.5 T), Uq{T) = Uo{0){l-T/Tef 
for core pinning, and to correct the value Uq (0) = 4300 K 
obtained there for the larger density of columns in those 
cases where B^ > 0.5 T. The difference between the low- 
field data and these values, represented by dashed lines 
in the Fig. |l|(b), again corresponds to temperature times 
the configurational entropy. For B^ — 1 and 2 T, it is 
found that TS{B = /SB^/e) ~ T, indicating that La ~ s. 
The logarithmic field dependence at T = 74 and 78 K, 
as well as the magnitude of the logarithmic slopes also 
agree with Lq/s « 1 — 2. However, since the logarith- 
mic slopes should, in the representation of Fig. |l|(b)), be 
equal to T, the predicted temperature dependence is not 
well obeyed. In fact, for the data at T = 70 K it is dif- 
ficult to identify any substantial field range over which 
TS ex ln{l/B). This is in contrast to data for B^ = 0.5 T 
where a low-field logarithmic dependence can be identi- 
fied at all investigated temperatures .liZlll3 Hence, accurate 
values for Hint can be determined only for small B^ (see 
Fig. H for Brj, = 0.5 T). At low B^, the entropy contribu- 
tion to the free energy is unimportant. The interaction 
field then delimits the low field regime where vortices are 
pinned independently on the columnar tracks, from the 
high field regime where intervortex repulsion entails the 
presence of unpinned ("free" or "interstitial") vortices. 
Results for larger matching field can only be obtained in 
a restricted temperature range Tirr < T < Tc{l — Gi) 
(see Fig. || for B^ — I T). Here, the interaction field sep- 
arates the regime of individual vortex pinning from high 
fields at which the intervortex repulsion determines the 
optimum pinned configuration but does not necessarily 
lead to the presence of "free" vortices. 

We note that, in the analysis presented in Fig. E|(b), the 
assumption of Uqs ss eo(0)s(l— T/Tc) for electromagnetic 
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FIG. 4. In-phase fundamental transmittivity T^ (a) and 
third harmonic amphtude \Th3\ (b) measured on an optimaUy 
doped Bi2Sr2CaCu208+i single crystal {Tc = 90.0 K), irra- 
diated with 5.8 GeV Pb ions to a matching field B^ = 0.5 
T. The measurement frequency was 7.753 Hz and the ac field 
amplitude hac = 6 Oe. DC field values are as indicated. 

pinning of vortices yields the rather unhkely result that 
the entropy contribution to the free energy rapidly de- 
creases as function of temperature. 



B. Relation with the phase diagram 

In order to correlate the vortex arrangement on the 
columnar tracks as revealed by the reversible magneti- 
zation measurements with the phase diagram, we have 
performed measurements of the ac transmittivity as de- 
scribed in section p^IB| . Typical results for the in-phase 
fundamental and third harmonic amplitude are shown in 
Fig. I for optimally doped BiaSraCaCuaOg+a {Tc = 90.0 
K) with B^ = 0.5 T£il The corresponding fields H„r{T) 
at which a third harmonic response is first detected are 
plotted in Fig. 0. Note that due to the different effec- 
tive measurement frequencies, / = 7.753 Hz in the ac 
shielding experiment versus / ^ 10^^ Hz for magneti- 
zation measurements using the SQUID, the iJirr~values 
plotted in the Figure lie much above the apparent Hirr- 
values extracted from the closing of the magnetic hystere- 
sis loops or merger of zero-field-cooled and fifild-pcooled 
data obtained with a SQUID magnetometer .eIEtEI The 
same is indicated in Fig. ^(a), which shows the position 
of Hirr determined from transmittivity measurements on 
the magnetization curve. The magnetization data show 
a break at Hirr , reflecting the inadequacy of moving 
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FIG. 5. Irreversibility field Hirr{T), at which a third har- 
monic ac response |Th3| can first be measured upon cooling 
the crystal (the same as in Fig. Wl) in a magnetic field (o). The 
roman numerals indicate the field regimes described in the 
text. Also indicated are the fields Hint(T) determined from 
reversible magnetization measurements on the same crystal 

sample SQUID magnetometry for measurements in the 
regime Ha ^ Hirr- Another reason for the much lower 
Hirr (T) values of Ref. p^|p^ is the fact that those authors 
used lightly underdoped crystals. While T,, can be similar 
for underdoped and overdoped samples, it was recently 
shown that the irreversibility lines for the irradictted ma- 
terial vary widely as function of oxygen content. l3 

The irreversibility field diplays three distinct regimes 
as function of temperature. Upon lowering tempei;a|turej 
from Tc, Hirr{T) first increases exponentially (I);EH3&e£| 
at r ~ 77 K, Hirr increases sharply (II), before bending 
over to an approximately linear temperature dependence 
below 70 K (HI). The same behavior is found for other 
matching fields, see e.g. Fig. pi for B^ = 1 T; for 
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FIG. 6. Irreversibility field pffirr(T) determined for lightly 
overdoped Bi2Sr2CaCu2 08+«,E3 with B^ = 1T (o). The crys- 
tal was chosen so as to have Tc = 83.1 K comparable to the 
sample used to obtain the JPR peak data of Ref. hll (also 
shown (o)), as well as to the crystal of Ref. Eq. The charac- 
teristic fields Bi (■) at which the c-axis conductivity starts to 
increase, and Bh ( D ) at which it starts its final decreaseEj 
are displayed as filled and open squares, respectively. The 
interaction fields Hi„t{T) determined from reversible magne- 
tization measurements are drawn as filled circles (•). 
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FIG. 7. Dependence of Hi,.r{T) on matching field. Dis- 
played are the values fj,(iHirr{T)/ B^ for a series of crystals 
grown at the University of Tokyo, irradiated with 5.8 GeV 
Pb ions to matching fields between 0.02 T and 4 T; data are 
plotted as function of T /Tc, where Tc is the value after irra- 
diation {e.g. Tc = 86.2 K for B^ = 2 T). 

comparison, the same Figure shows characteristic fields 
extracted from c-axis conductivity and JPR measure- 
ments. It is seen that the different regimes of the irre- 
versibility line Hirr(T) ate, faithfully reproduced by the 
position of the JPR peak;tll for suitably chosen measure- 
ment frequency (here, 41 GHz) the two lines lie very close 
together, suggesting that the onset of linear vortex diffu- 
sion (and hence, the disappearance of long-range a5-plane 
phase coherence) indicated by Hirr{T) is intimately re- 
lated to the drop in c-axis phase coherence expressed by 
the JPR peak. 

Also plotted in either Figure are the TJ^t-data for the 
same crystals. For the optimally doped crystal {B^ = 0.5 
T) of Fig H, it was possible to obtain i?i„t-values di- 
rectly from the reversible magnetization for fields below 
the irreversibility field by exploiting the torque data of 
Ref. ^. The reversible magnetization curves in this case 
were measured with the magnetic field applied at a sub- 
stantial angle with respect to the c-axis and the column 
direction. It was shown previously that in such a con- 
figuration, the reversible magnetization depends only on 
the field component H^ parallel to the c-axis; within ex- 
perimental error, the magnetization curve MreviHz) is 
the same regardless of the direction of the total magnetic 
field. Simultaneously, applying the field at an angle low- 
ers the irreversibility field Hirr to a much lower value 
than when the field is aligned with the columnar defects.El 

It appears that the field Hint intercepts the irreversibil- 
ity line exactly at the point fioHa ~ g^"^ where Hirr 
starts its steep rise (Fig g). This strongly suggests that 
the reduction of vortex mobility causing the increase in 
Hirr (T) is the consequence of intraplane vortex repulsion 
and the reduction of the number of sites available to a 
given vortex. The Hint data nicely coincide with the re- 
coupling field obtained in Refs. ^ 12. For B^ = 1 and 2 
T, reliable 7Ji„f-data could not be obtained in the tem- 
perature interval spanning the irreversibility line. The 



data in Fig. g show that for temperatures well in excess 
of Hirr, Hint dccrcascs with increasing T. The temper- 
ature dependence of the interaction field, reminiscent of 
that of the high-field part of the irreversibility line, seems 
to extrapolate to the field B^ where the c-axis conduc- 
tivity starts its ultimate decrease. Ej Note that Hint never 
increases beyond the field fioHa ~ ^B^ at which the re- 
coupling transition is expected. This particular fraction 
of the matching field, ^B^, seems to be quite robust. 
Figure |^, which displays the evolution of HirriT) with 
Bfj,, shows that for 0.1 T < B^j, < 4 T the steep rise of 
the irreversibility line always starts near //Or^a ~ g-B^- 
This is in contradiction tO|-^umerical worktZI and pre- 
vious magnetization studiesta which reported significant 
features in the phase diagram at n^Ha ~ hB^p. Fig- 
ure shows that, whereas for the values of the matching 
field Bcf, — 1 aa d 2 J T, most commonly investigated in 
the literature, Ilirt^EZl the fraction -B^/S coincides with 
the crossover from regime (II) to regime (III) in Hirr{T), 
this is by no means true for other values of B^. The 
results also show that the reported scaling of the irre- 
versibility line with B^ reported in Ref. ra is by no means 
to be taken strictly. 



IV. DISCUSSION 

A number of significant results emerge from the above. 
The most notable is that the introduction of amorphous 
columnar defects into the matrix of the layered supercon- 
ductor decreases the free energy in the London regime 
not only through vortex pinning, but also through an 
additional configurational entropy, which appears when 
a sufficiently large number of extra sites become available 
to the pancake vortices. This extra contribution can be 
comparable to or larger than the pinning energy. Pre- 
viously presented values of the pinning energy extracted 
from the reversible magnetiZ|ation of crystals with rela- 
tively high matching ficldstJ Il3£3 should thus be consid- 
ered with reserve. More accurate measurements of the 
pinning energy for crystals with lower B^ were recently 
presented in Ref. Qq. The pinning energies obtained in 
this paper are in agreement with those of Ref. |lj apd with 
the vortex-core pinning model, Uq ~ eoico/'^S,)'^^ This 
means that even in the vortex liquid state, vortices are, 
certainly in the regime B <C i?0, strongly bound to the 
columnar defects. Thermally induced vortex wandering 
from the tracks would have resulted in a pinning energy 
contribution Uq/^o that is expcpentially small with re- 
spect to the total magnetization.E3 The neglect of the en- 
tropy contribution also accounts for the discrepancy be- 
tween experimental results and the numerical simulations 
of Ref. pl^. The curves calculated in Ref. O are close to 
what would be expected in the (near-) absence of an en- 
tropy contribution (see for example data at i?^ = 0.2 T). 
Conversely, the use of Eq. (^ overestimates the entropy 
contribution in many cases and yields pinning energies 
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FIG. 8. Entropy contribution to the free energy of 
heavy-ion irradiated hghtly overdoped Bi2Sr2CaCu2 08+5 
(Tc = 83 K) as function of matching field (o), for T = 74 
K and applied field jJ,oHa =0.1 T. The entropy rises steeply 
between B^ = 0.5 and 1 T, from where it increases loga- 
rithmically as function of -B^. This behavior is in qualitative 
agreement with the dissociation of individual vortex lines into 
stacks of length La at low defect density (few available sites 
per vortex), for which the prediction TS oc B^lnBcj, is indi- 
cated by the curved full line; at higher matching fields, the 
length of the typical stack has decreased to the layer spacing 
s — the dependence then becomes TS oc In _B^ corresponding 
to the configurational entropy associated with the possibility 
for each pancake vortex to occupy difi^erent columns (dashed 
straight line). The corresponding values of La are indicated 
by the full triangles (the dotted line is a guide to the eye)n 
The parameter values used are eqs = 1500 K, Bcr ~ 0.4 tJ^ 
i.e. 7 = 55.Ej 

that are manifestly too smalla 

The measurement of the reversible magnetization as 
function of column density allows a direct determination 
of the extra entropic contribution TS (Fig. g), and its 
dependence on field and temperature. The mechanisms 
by which the free energy is lowered becomes apparent if 
one considers the rather peculiar dependence of TS on- 
the density of sites at constant field, illustrated in Fig. ||. 
At low B^, TS remains small in comparison with other 
terms contributing to the free energy. This can also be 
inferred directly from the near-equality of the low- and 
high field limits of the logarithmic-slope of the reversible 
magnetization for B^ < 0.5 T.EZHla At high matching 
field, TS follows a logarithmic dependence on ion dose. 
We interpret this behavior as follows. At low irradiation 
doses, the increase in entropy with increasing columnar 
defect density arises from the dissociation of vortex lines 
into shorter and shorter stacks of length La localized on 
different columns; the entropic part of the free energy (in 
the vortex liquid phase) then obeys Eq. (0) , with /? ~ i . 
As a result, TS increases as B^lnB^, which explains its 
rapid rise between B^ = 0.5 T and 1 T and the quite im- 
portant change in the low-field magnetization between 
these matching fields. Once the length of a typical stack 
has decreased to a single pancake vortex, this mechanism 
can no longer be operative. The entropy now increases 
as InS^, corresponding to the possibility for each single 



pancake to occupy different columnar defects. The evo- 
lution of the entropy with increasing matching field can 
be described by the optimisation of the total free energy 
gain, obtained from the balance between the entropy gain 
and the loss in intervortex interaction energy. At small 
-00, there are few near lying columnar defect sites per 
vortex line — the access to farlying sites is prohibited 
by intervortex repulsion. Hence, vortices are on average 
well-aligned on the same site. For increasing B^, more 
and more nearlying defect sites become available, and 
the free energy can be lowered through the dissociation 
of vortex lines, first into stacks of pancakes located on 
adjacent columns, and finally into single pancakes wan- 
dering between available column sites. 

The same explains the field dependence of TS (Figs.|l| 
and H). As more and more vortices are added, fewer 
sites are available per vortex. First, the entropy contri- 
bution drops (logarithmically) because fewer defect sites 
are available to each pancake. At higher fields, pancakes 
belonging to the same vortex line become aligned on the 
same site, be this on a columnar defect or on an intersti- 
tial site. This leads to a reduction in vortex mobility, and 
to tlie, mcrease of the c-axis critical curr.eiit measured by 
JPf^y^EJ and of the c-axis conductivity.ElO'El 

It is to be stressed that in the presence of columns the 
low-field entropy contribution to the free energy is much 
larger than that in unirradiated crystals, at odds with 
the prediction of Eqs. (y,|^).C3 A qualitative explanation 
for this is that, whereas in unirradiated crystals one has 
one local energy minimum per vortex lattice unit cell, 
the presence of a columnar defect-density nj, > B/^q 
amounts to the presence of several energy minima per 
cell, permitting the "roughening" of vortex lines {i.e. 
their spreading over different columnar defects). Entropy 
decreases as function of B until at i? ~ ^B^ it attains a 
value close to that measured in unirradiated samples, cor- 
responding to approximately one site per vortex. From 
this field onwards, it becomes difficult to obtain a reli- 
able estimate of the alignment length La from reversible 
magnetization. This is because of the "error bar" intro- 
duced by the large valu e of the Ginzburg number Gi. 
We have seen in Section [III A 2 there is an uncertainty 
of the order O.lksT/^os when one determines the rel- 
ative importance of the different terms contributing to 
dMrev/dhiB. In particular, one cannot estimate the 
entropy contribution to better than O.lksT/^QS, corre- 
sponding to a lower bound La = 10s. The same holds 
for unirradiated crystals. Whereas, in the irradiated su- 
perconductor, the dissociation of vortex lines into stacks 
of pancakes (within a unit cell) certainly lowers inter- 
plane correlations of the phase of the superconducting 
order parameter, it does not necessarily imply the loss 
of long-range phase coherence. This is likely destroyed 
near the irreversibility line, as linear activated motion of 
(pancake) vortices becomes possible. 

The questions of pancake vortex alignment, phase 
coherence, and the recoupling transition reported in 
Refs. pT|-p^ are intimately related. From reversible mag- 
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netization, we have determined the field Hi„t at which 
intervortex repulsion determines the vortex pancake ar- 
rangement in the presence of columnar defects. This is in 
very good agreement with the recoupling field obtained in 
Refs. 



11 



. 12, and coincides with the rapid rise of Hirr {T) . 
It seems that the same intraplane vortex repulsion re- 
sponsible for vortex alignment also decreases the vortex 
mobility, expressed by the peak in the critical current 
density measured in Refs. |l3 and q1|, and by the increase 
in Hirr{T) into regime (II) at \B^. As for the increase in 
vortex mobility and the decrease of c-axis conductivity 
at the crossover of Hirr from regime (II) to the high-field 
region (III), it is likely that it corresponds to the field at 
which the number of free vortices becomes comparable to, 
or exceeds the number of vortices trapped on a colum- 
nar defect. This field can be roughly estimated as that 
at which Up + TS (Figs. |l|(b) and y(b)) has decreased 
to half its value at Hint- For i?^ > 1 T, it amounts to 
roughly the value hB^f, cited in Ref. |l3|; for lower match- 
ing fields it lies closer to B^. In all cases, it lies close to 
the crossover in question. Thus, we identify regime (III) 
of Hirr{T) as that in which the number of free vortices 
exceeds the number of trapped ones, and in which vortex 
mobility is determined by plastic vortex creep .E2 



V. CONCLUSION 

Reversible magnetization measurements have been an- 
alyzed in order to determine the relative importance of 
pinning energy and entropy contributions to the free 
energy of heavy-ion irradiated Bi2Sr2CaCu208+5 single 
crystals. We have found that in unirradiated crystals the 
entropy contribution to the free energy in the London 
regime is relatively minor. This is also true in irradiated 
crystals with small matching fields B^ < 0.5 T, but not 
for irradiated crystals with large B^ > 1 T and fields 
B <C B^. Then, the configurational entropy contribu- 
tion to the free energy is greatly enhanced with respect 
to its value in unirradiated crystals, and can easily ex- 
ceed the pinning energy of a vortex on a columnar defect 
by more than a factor 2. The large entropy contribu- 
tion at high i?0 is a consequence of the dissociation of 
vortex lines into individual pancake vortices localized on 
different columnar defect sites. As B/B^ is increased, in- 
tervortex repulsion causes the pancakes belonging to the 
same stack to align onto the same columnar defect. It was 
possible to obtain a lower bound La ~ 10s on the length 
over which pancakes are aligned at _B « i?^. Vortex in- 
teractions become important at the field Hint ~ hB,^,', at 
this field, which delimits the regime of pinning of indi- 
vidual vortices, and which corresponds to the recoupling 
field measured by Josephson Plasma Resonance, the irre- 
versibility field sharply increases. Above Hint, the repul- 
sive interaction between vortices causes both the vortex 
mobility to decrease and pancake alignement to increase. 
At higher fields ^ |-B^, free vortices outnumber those 



that are trapped on a columnar defect. This causes the 
decrease of c-axis resistivity and a second crossover of the 
irreversibility field, to a regime where it is determined by 
plastic creep. 
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